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摘 要 : 当前 基于 中 等 空间 分 辩 率 时 序数 据 的 农作物 种 植 结构 提取 成 为 研究 热点 ,但 农作物 季 相 节律 特征 在 不 同 
气候 背景 下 存在 较 大 差异 ,绿洲 作为 干旱 区 具有 明显 小 气候 效应 的 生态 景观 ,其 农作物 种 植 结构 的 还 感 提取 具有 


较 强 的 典型 性 和 代表 性 。 选 取 宁 夏河 套 平原 绿洲 典型 区 域 , 通 过 构建 高 分 


号 (GF -1AWFV ) 时 间 序 列 数据 ,结合 


不 同 作物 耕作 方式 及 生长 物候 ,分 析 不 同 作物 在 整个 生长 季 内 的 归 一 化 植被 指数 (NDVI) 和 归 一 化 水 体 指数 ( ND- 


WI) 的 时 间 序 列 特征 ,构建 不 同 决策 树 提取 研究 区 农作物 种 植 结 构 信 息 , 并 验证 了 不 同方 法 的 适用 性 。 


结果 表明 ， 


对 具有 明显 小 气候 效应 的 干旱 区 绿洲 ,利用 时 间 分 辨 率 和 空间 分 辩 率 都 较 优 的 GF1-WFV 时 间 序 列 数据 ,对 其 农 作 


物种 植 结 构 进 行 遥 感 提取 具有 较 强 的 实用 性 和 代表 性 。 
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农作物 种 植 结构 是 指 一 个 地 区 内 主要 农作物 的 
类 型 和 分 布 。 和 掌握 作物 的 种 植 结构 情况 ,可 为 农业 
精准 管理 提供 重要 的 参考 依据 。 遥 感 技术 具有 宏观 
性 、 综 合 性 和 动态 性 的 特点 ,是 获取 区 域 农作物 布局 
和 组 成 的 重要 途径 "“。 传 统 的 利用 遥感 技术 提取 
农作物 种 植 结构 多 基于 单一 各 感 影 像 进 行 提取 ,该 
方法 操作 简单 ,但 往往 难以 确定 种 植 结 构 “ 最 佳 识 
别 期 "的 遥感 影像 中 ;基于 时 间 序 列 数据 的 作物 类 
型 提取 起 初 多 见于 以 MODIS 为 代表 的 低空 间 分 辩 
BDO ,其 在 大 区 域 尺度 上 农作物 种 植 结构 的 
提取 取得 了 较 好 的 效果 ;在 中 等 分 辨 率 数 据 源 方面 ， 
以 Landsat 卫星 为 代表 , 受 时 间 分 辨 率 的 制约 ,多 集 
中 在 基于 几 个 时 相 的 作物 种 植 结构 提取 研究 '" ; 随 
着 遥感 数据 源 的 不 断 丰 富 ,利用 高 空间 或 高 时 间 分 
辨 率 遥 感 影 像 提取 农作物 种 植 结 构 已 成 为 国内 外 研 
究 热点 。2013 年 4 月 26 日 我 国 高 分 一 号 卫星 
(GF - 1) 发射 运行 ,通过 4 台 宽 幅 相 机 视 场 拼接 ,在 
16 m 的 空间 分 辨 率 上 ,实现 了 2 d 的 重 访 周期 ,大 大 
提高 了 中 等 空间 遥感 数据 的 时 间 分 辩 率 。 

不 同 作物 的 耕作 方式 不 同 , 在 发 育 过 程 中 时 间 
和 生物 量 上 存在 一 定 差 异 , 即 物候 特征 不 同 ” , 单 
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一 时 相 遥 感 影像 数据 中 不 同类 型 农作物 光谱 发 生 交 
又 或 重 释 ,或 农作物 物候 差异 不 明显 ,因而 使 得 农 作 
物种 植 结构 提取 发 生 混淆 和 模糊 ,基于 多 时 相 
遥感 影像 数据 构建 时 间 序 列 已 成 为 作物 分 类 研究 的 
热点 。 基 于 野外 实地 调查 数据 构建 作物 特征 指数 时 
间 序 列 ,结合 作物 物候 特征 引入 决策 树 算法 ,实现 作 
物种 植 结构 的 半自动 化 提取 ,是 当前 大 面积 作物 分 
布 信息 监测 研究 的 重要 方向 。 决 策 树 算 法 是 
一 种 非 参 数 化 , 非 线 性 的 监督 分 类 法 ,具有 构造 简单 
灵活 ,精度 较 高 和 结果 语义 明晰 .可 解释 性 强 等 特 
点 ,广泛 应 用 于 土地 科学 和 地 学 研究 各 领域 ,尤其 是 
在 遥感 图 像 分 类 .土地 评价 .土地 利用 类 型 的 提取 等 
Tria" 。 目 前 ,比较 成 熟 的 决策 树 构 建 方法 有 
QUEST ID3 、C5. 0, CART 和 曲线 特征 等 ~ ,其 
中 ,基于 曲线 特征 构建 决策 树 和 CART 决策 树 是 两 
种 常用 的 决策 树 分 类 方法 。 刘 晓 娜 等 "基于 地 物 
的 不 同 特征 曲线 构建 决策 树 , 实 现 了 对 橡胶 林 的 遥 
感 识别 。 常 布 辉 等 2 光 利 用 曲线 振幅 和 相位 差 构建 
决策 树 ,并 结合 监督 分 类 对 河套 灌区 沈 乌 灌 域 耕地 
的 提取 。 于 文 婧 等 2 建立 研究 区 内 典型 地 物 的 
NDVI 时 间 序 列 曲线 ,然后 对 土壤 信息 丰富 的 3 月 多 
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光谱 影像 进行 主 成 分 变换 ,选取 第 1 主 成 分 (PC1 ) 
作为 光谱 特征 参数 ,最 后 建立 CART 决策 树 实现 了 
对 平 罗 县 土地 利用 的 遥感 分 类 。 这 些 研 究 均 表明 ， 
决策 树 方 法 分 类 结果 清晰 准确 , 旦 精度 高 于 其 他 分 
类 方法 。 

我 国 是 世界 上 绿洲 面积 最 大 、 类 型 最 丰富 .人 研究 
也 最 多 的 国家 ,绿洲 在 干旱 区 社会 经 济 发 展 过 程 中 
起 着 至 关 重 要 的 作用 , 且 作 为 干旱 区 具有 明显 小 气 
候 效应 的 生态 景观 “1 ,其 农作物 种 植 结构 的 遥感 
提取 具有 较 强 的 典型 性 和 代表 性 。 以 地 处 宁夏 河套 
绿洲 中 段 的 贺 兰 县 和 平 罗 县 为 研究 区 ,对 GFI - 
WFV 影像 数据 进行 预 处 理 , 通 过 波段 运算 得 到 植被 
指数 (NDVI) 、 水 体 指数 (NDWI) 的 序列 数据 ,计算 
和 提取 反映 该 区 主要 作物 物候 特征 的 时 序 参 数 , 分 
析 作 物 时 序 特征 与 耕作 信息 的 关联 ,分 别 采用 基于 
时 间 序 列 曲线 特征 的 决策 树 分 类 方法 和 CART 决策 
树 分 类 方法 提取 该 地 区 的 农作物 种 植 结 构 , 并 分 析 
对 比 了 不 同 决策 树 的 优 缺 点 和 适用 性 。 


1 研究 区 概况 


研究 区 位 于 贺兰山 与 鄂尔多斯 高 原 之 间 , 东 临 
黄河 , 西 靠 贺 兰 山 ,地 理 坐 标 为 106°*48”~ 106°16'E， 
38°26" ~39?°00'N, 是 在 新 生 代 断 陷 盆 地 基础 上 发 育 
的 堆积 平原 2 , 自 西向 东 由 山 前 洪 积 倾斜 平原 、 冲 
洪 积 平原 和 河 湖 积 平原 组 成 ,地 势 平缓 低洼 。 该 区 
属于 温带 大 陆 性 干旱 气候 , 春 旱 多 风 , 夏 热 多 雨 , 秋 
季 短 暂 ,冬季 干 冷 ,年 平均 气温 8.8 C ,年 平均 降水 
Æ 173 mm ,降水 主要 集中 在 7 一 9 月 ,年 平均 蒸发 量 
1755 mm, ARBRE WES 月 。 研 究 区 日 照 充 足 ， 
早晚 温差 较 大 ,总 热量 丰富 ,每 年 无 霜 期 较 长 ,在 适 
合作 物 和 后 长 的 4 一 9 H S10 % 的 积温 在 3 200 ~3 400 
C ,能 满足 小 麦 , 水 稻 等 喜 温 作物 的 生长 需要 。 并 且 
210 C 的 积温 日 期 也 与 无 霜 期 十 分 吻合 ,太阳 辐射 
较 高 ,多 年 平均 日 照 时 间 为 2 800 ~3 100 h ,无 霜 期 
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图 1 宁夏 河套 绿洲 作物 种 植 结构 实测 数据 示意 图 
Fig.1 Measured data of crop planting structure in the 


Hetao Oasis in Ningxia 


为 164 d 左右 ,这 些 条 件 均 有 利于 作物 的 正常 生长 。 

宁夏 河套 平原 绿洲 作为 中 国 西北 地 区 传统 的 重 
要 粮食 生产 基地 ,长 期 以 来 形成 了 以 粮食 作物 为 主 
的 种 植 结构 体系 (图 1)。 主 要 作物 有 水 稻 (4 一 9 
月 ) 玉米 (4 一 9 月 )、 小 麦 (3 一 7 A) Aia (3—11 
月 ) 等 ,其 生育 期 见 表 1。 河 套 绿 洲 是 依靠 黄河 灌溉 
而 成 的 外 流 型 绿洲 ,水利 开 发 对 河套 绿洲 农业 具有 
决定 性 意义 ,宁夏 河套 平原 早 在 2000 多 年 前 的 秦 、 
汉 时 代 就 引用 黄河 水 进行 农田 灌溉 ,如 今 有 14 条 主 
灌溉 汇 和 多 条 排水 沟 , 灌 排 水 系统 发 达 。 宁 夏河 套 
绿洲 由 于 地 处 我 国 干旱 地 区 , 莹 发 量 是 降水 量 的 9 
倍 之 多 ,具有 明显 的 小 气候 效应 ,又 是 典型 的 灌溉 农 
业 区 ,所 以 其 农作物 类 型 的 遥感 提取 具有 较 强 的 典 
型 性 和 代表 性 。 


2 数据 来 源 与 处 理 


2.1 遥感 数据 获取 与 处 理 
高 分 一 号 卫星 于 2013 年 4 月 26 日 发 射 , 其 中 


表 1 宁夏 河套 绿洲 主要 作物 生育 期 


Tab.1 Phenological period of the main crops in Ningxia Hetao oasis 


主要 3 月 4 月 5H 6 H 7A 8 A 9 月 10 月 
农作物 上 中 下 上 中 下 上 中 F 上 中 F 上 中 下 上 中 下 上 中 F 上 
KAR 灌水 灌水 DEE oye 拔节 拔节 孕穗 抽穗 抽穗 灌浆 灌浆 成 熟 成 熟 
玉米 出 苗 出 苗 拔节 灌浆 灌浆 FLA FLA 成 熟 成 熟 

小 麦 出 苗 aT BE 拔节 拔节 抽穗 抽穗 灌浆 灌浆 成 熟 成 熟 


EEG Fa) XU Fl 


FX Fl XR XU Fal X XE 


刘 雅 清 等 :基于 GF -1AWFV 时 间 序 列 的 绿洲 作物 类 型 提取 


表 2 时间 序列 数据 信息 
Tab.2 Information of the time series data 
序号 传感器 类 型 和 产品 序号 影像 获取 时 间 | 序号 传感器 类 型 和 产品 序号 影像 获取 时 间 
1 GF1_WFV4_E106. 7_N38.5_L1A0002263952 2017 -03 -26 8 GF1_WFV2_E107.4_N39.3_L1A0002473378 2017 -07 -10 
2 GF1_WFV4_E107.3_N38.3_L1 A0002331148 2017 -04 -28 9 GF1_WFV3_E105. 8_N38.9_L1A0002525744 2017 -08 -04 
3 GF1_WFV3_E105. 9_N38.9_L1 A0002344768 2017 -05 -06 10 GF1_WFV3_E106. 1_N38.9_L1A0002545275 2017 -08 -16 
4 GF1_WFV3_E106. 2_N38.9_L1A0002379084 2017 -05 -26 11 GF1_WFV4_E106. 6_N38.5_L1A0002580944 2017 -09 -02 
5 GF1_WFV4_E106. 7_N38.5_L1A0002424156 2017 -06 - 16 12 GF1_WFV2_E105. 9_N39.3_L1A0002649081 2017 - 10 -04 
6 GF1_WFV3_E106. 6_N38.9_L1A0002449838 2017 -06 -28 13 GF1_WFV2_E107. 5_N39. 3_L1A0002671474 2017-10-12 
7 GF1_WFV3_E106. 1_N38.9_L1A0002464407 2017 -07 -06 14 GF1_WFV2_E107. 1_N39.3_L1A0002755542 2017 -11 -10 


搭载 了 4 台 WFV 多 光谱 相机 ,4 台 相 机 组 合 扫描 幅 
宽 800 km , 星 下 点 的 分 辩 率 16 m, WEV 传感器 共 设 
置 4 个 波段 ,光谱 范围 为 0.45 ~0.89 pm( 蓝 光 
0.45 ~0.52 um , 绿 光 0.52 ~0.59 hm, 红 光 0.63 ~ 
0.69 um, 近 红外 0.77 ~ 0. 89 um) , 重 访 周期 为 
2 d。 本 研究 自 2017 年 3 一 11 月 共 获 取 14 景 影像 ， 
履 盖 研究 区 所 有 作物 的 生长 周期 ( 表 2) 。 选 取 影 像 
时 ,确保 研究 区 影像 云 量 最 小 ,影像 的 时 间 尽 量 靠近 
作物 特殊 物候 期 。 利 用 ENVI 5. 3 遥感 影像 处 理 软 
件 对 多 时 相 的 GF - 1/WEV 影像 进行 辐射 定 标 ,得 
出 大 气 上 行 辐射 亮度 值 ,并 运用 FLAASH 模型 进行 
大 气 校正 ,所 需 的 基本 参数 数据 来 源 于 影像 头 文件 
和 中 国资 源 卫 星 中 心 。 通 过 研究 区 矢量 边界 对 多 时 
相 影 像 进行 拼接 和 裁剪 ,得 到 人 研究 区 影像 集 。 最 后 
根据 野外 实测 GPS 控制 点 ,采用 二 次 多 项 式 模型 ， 
对 影像 数据 进行 几何 精 校正 ,校正 误差 控制 在 0.5 
个 像 元 以 内 。 
2.2 野外 采样 数据 获取 与 处 理 

野外 调研 时 发 现 ,研究 区 内 地 物 主要 有 湖泊 鱼 
塘 、 林 地 建筑 用 地 和 荒地 天然 草 地 水稻、 玉米 小 
麦 . 首 蒂 和 和 葡萄、 枸杞 等 少量 其 他 作物 。 根 据 研究 区 
内 作物 种 植 特点 ,主要 提取 水 稳 \ 玉 米 .小麦 和 首 蒂 ， 
其 他 农作物 分 布 面积 较 小 ,本 研究 忽略 。 据 此 ,将 研 
究 区 地 物 覆 六 类 别 分 为 水 稳 、 玉 米 .小麦 . 首 蒂 、 其 他 
植被 和 非 植 被 。 

利用 手持 CPS 对 主要 作物 进行 定位 调查 , 共 标 


到 NDVI 和 NDWI 时 间 序 列 通 感 数据 集 ,通过 训练 
样本 点 ,提取 不 同 作物 的 NDVI 和 NDWI 值 ,构建 
NDVI 和 NDWI 时 间 序 列 曲线 ;(3 分 析 不 同 植被 的 
肥 感 指数 时 间 序 列 数 据 特征 ,基于 时 间 序 列 曲 线 的 
振幅 和 相位 确定 作物 的 提取 阔 值 ,构建 决策 树 ,提取 
作物 种 植 结构 ;@ 将 训练 样本 作为 目标 变量 导入 
Clementine 软件 ,确立 水 稳 、 玉 米 和 小 麦 的 精确 分 类 
BO {EL PAE CART 决策 树 提 取 作 物种 植 结构 ;G5) 验 
证 对 比分 类 精度 ,分 析 讨 论 。 
3.2 指数 选择 与 计算 

参考 提取 农作物 种 植 结构 的 相关 文献 -和 
研究 区 主要 作物 的 光谱 特征 及 物候 特征 ,选取 归 一 
化 植被 指数 (NDVI) 和 水 体 指 数 (NDWI) 作 为 提取 
研究 区 农作物 种 植 结构 的 变量 参数 并 构建 时 间 序 
列 。NDVI 也 称 为 生物 量 指标 变化 ,可 使 植被 从 水 
和 土 中 分 离 出 来 ,这 在 一 定 程度 上 消除 太阳 高 度 角 、 
卫星 观测 角 、 地 形 、 云 以 及 大 气 条 件 的 影响 。 研 究 区 
水 域 范围 广 ,水 稻 种 植 面积 大 , 且 水 稻 和 玉米 的 生长 
周期 较为 相似 ,而 NDWI 可 以 突出 水 体 信 息 , 由 于 水 
体 的 反射 从 可 见 光 到 中 红外 波段 逐渐 降低 ,到 近 红 
外 范围 内 几乎 无 反射 ,而 植被 一 般 在 近 红 外 波段 的 
反射 率 最 强 , 因 此 NDWI 对 于 水 稻 区 分 十 分 有 利 。 
各 指数 计算 公式 见 表 3。 


表 3 特征 指数 计算 公式 


Tab.3 Formulas for calculating the characteristic exponents 


记 727 个 地 块 ,其 中 ,水 称 449 个 ,玉米 170 个 ,小 麦 
83 个 , 首 蒂 25 个 ,分 为 训练 样本 和 验证 样本 。 


3 研究 方法 


3.1 技术 路 线 
D 进行 实地 调查 ,定位 不 同 作 物 的 坐标 ;@) 对 
完成 预 处 理 的 GF - 1/WEV 影像 进行 波段 运算 ,得 


特征 指数 计算 公式 
NDVI NDVI =pyin ~ Pr/Pnir + PR 
NDWI NDWI=p¢ -pvr /pe +PNIR 


注 :pNm 为 近 红 外 波段 的 反射 率 ;p 为 红 光 波段 的 反射 率 ;pc 为 绿 光 
波段 的 反射 率 。 


3.3 2 种 特征 指数 时 序 构建 
3.3.1 时 间 序 列 曲线 构建 ”训练 样本 的 选择 对 分 
类 信息 的 获取 和 分 类 规则 的 建立 起 关键 作用 ,在 选 


THR 


WIRE E PES A FILAN: D 训练 样 
本 点 均匀 分 布 在 研究 区 范围 内 ,确保 选取 的 训练 样 
本 点 有 代表 性 ;@ 训练 样本 点 所 在 范围 地 块 面积 
大 ,作物 长 势 均 匀 ;(3) 训练 样本 个 数 占 总 样本 的 2% 
以 上 5 ,满足 完备 性 和 代表 性 两 个 基本 原则 。 基 于 
训练 样本 获取 不 同 作物 NDVI 和 NDWI 并 计算 其 均 
值 ,绘制 不 同 作物 NDVI 和 NDWI 时 序 特征 曲线 (图 
2 和 图 3)。 

3.3.2 作物 时 间 序 列 曲线 分 析 ”从 图 2 和 图 3 可 
以 看 出 ,研究 区 水 稳 在 灌水 期 即 4 月 下 旬 至 5 月 下 
旬 NDYVI 值 小 于 玉米 和 小 麦 ,期 间 其 NDWI 值 远大 
于 其 他 作物 ,5 月 下 旬 开 始 直 播 稻 出 苗 ,插秧 稻 移 
A, BO HANTERTH, E NDVI 值 增长 较 快 ,7 
月 进入 孕穗 抽穗 期 ,8 月 初 进入 灌浆 乳 熟 ,NDVI 值 
增 至 最 大 ,9 月 初 进入 灌浆 后 期 ,NDVI 值 逐 渐 下 降 。 
玉米 于 4 月 下 旬 播 种 ,5 月 中 旬 至 6 月 中 旬 进 入 拔 
节 期 与 开花 授粉 期 ,其 NDVI 值 增 速 较 快 ,在 7 月 初 
至 8 月 初 NDVI 值 较为 稳定 ,8 HRF RA FLA 
期 ,NDVI 值 逐渐 降低 。 玉 米 与 水 稳 的 NDVI 曲线 变 
化 趋势 较为 相似 。 小 麦 在 3 月 出 苗 ,4 月 进入 分 组 
期 ,在 5 月 进入 拔节 孕穗 期 , 在 此 期 间 小 麦 快 速生 


i: j 
5 0.8 L = EK 
© 0.7 上 一 小麦 
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图 2 各 地 类 样本 的 NDVI 时 间 序列 曲线 


Fig.2 Time series of NDVI of different sample plots 
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图 3 各 地 类 样本 的 NDWI 时 间 序 列 曲线 
Fig.3 Time series of NDWI of different sample plots 
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长 ,其 NDVI 值 远大 于 水 稻 和 玉米 ,6 A RAFA 
期 ,其 NDVI 值 迅 速 下 降 ,7 月 初 进入 收获 期 ,其 ND- 
VI 降 至 最 低 , 曲 线 出 现 波 谷 , 小 麦收 割 后 受 田间 杂 
草 生长 的 影响 ,其 NDVI 曲线 在 8 月 中 旬 会 出 现 第 
二 个 峰值 。6 月 中 旬 小 麦 灌 浆 期 消 小 麦 灌 浆 水 ,6 月 
下 旬 小 麦 消 麦 黄 水 ,期 间 小 麦田 块 的 NDWI 值 有 明 
显 的 上 升 趋势 。 首 蒂 返 青 时 间 早 ,是 年 内 会 多 次 重 
复出 现 生 长 一 刘 制 循环 , 首 蒂 田 块 现存 生物 量 也 因 
此 具有 逐渐 积累 一 高 峰 期 一 收获 后 的 最 低 时 期 的 变 
化 过 程 , 导 致 在 生长 季 内 NDVI 呈现 出 由 低 逐 渐 到 
高 .由 高 迅速 降低 的 循环 变化 特征 , 且 在 每 次 浊 割 后 
首 蒂 田 块 会 进行 灌水 施肥 的 农业 活动 ,所 以 其 ND- 
WI 也 呈现 出 由 低 逐 渐 到 高 .由 高 迅速 降低 的 循环 变 
化 特征 。 
3.4 构建 决策 树 分 类 规则 

从 土地 覆 被 角度 来 看 ,水 体 类 与 非 水 体 类 以 及 
植被 类 与 非 植被 类 ,特定 的 地 物 具 有 其 特定 的 时 序 
特征 ,只 要 依据 地 物 时 序 特 征 合 理 地 选择 地 物 特 征 
序列 组 合 , 就 能 把 目标 地 物 与 其 他 地 物 区 别 开 ™?。 

已 有 研究 表明 5 ,NDWI 大 于 0 即 为 水 体 ,ND- 
VI 大 于 0.3 为 植被 ,小 于 0.3 为 非 植被 。 研 究 区 河 
流 、 湖 泊 、 鱼 塘 等 水 域 面积 较 大 ,对 后 续 提取 水 稳 干 
扰 较 大 。 经 调查 ,研究 区 3 月 无 灌溉 活动 ,为 排除 干 
扰 , 可 利用 研究 区 3 月 26 日 影像 的 NDWI 值 提 取出 
水 体 。 研 究 区 6 月 中 旬 为 所 有 作物 生长 旺盛 期 ， 
NDVI 值 较 大 ,利用 6 月 16 日 影像 的 NDVI 值 可 以 
区 分 植被 与 非 植被 ,将 建筑 用 地 和 未 利用 荒地 等 提 

首 蒂 生长 季 内 NDVI 呈现 出 由 低 逐 渐 到 高 、 由 
高 迅速 降低 的 循环 变化 特征 。 鉴 于 此 ,可 采用 二 次 
差分 法 ”提取 极 值 计算 波峰 个 数 的 方法 提取 首 
董 。 由 于 二 次 差分 法 对 极 值 的 提取 非常 敏感 ,能 够 
提取 每 一 个 细小 的 极 值 ,为 了 避免 错误 极 值 对 提取 
波峰 产生 干扰 ,该 研究 在 判断 峰值 时 确立 波峰 NDVI 
值 大 于 0.5 进行 波峰 的 取舍 。 

水 稻 于 4 ARABS 月 下 旬 开 始 灌水 ,于 5 月 
23 日 水 稳 田 灌水 已 基本 结束 。 灌 水 后 水 稻田 的 
NDVI 值 大 于 其 他 作物 的 NDWI 值 ,可 基于 5 月 26 
日 研究 区 影像 NDWI 值 提取 研究 区 水 称 种 植 范围 ， 
对 水 稳 和 芦 全 的 区 分 选用 10 月 12 日 的 NDVI 值 。 

小 麦 在 3 月 播种 ,至 5 月 初 已 到 分 莫 至 拔节 期 ， 
而 此 时 直播 稻 和 玉米 刚刚 出 苗 , 移 栽 稻 还 未 移 栽 , 影 
像 NDVI 值 远 小 于 小 麦 ,小麦 在 6 月 下 旬 进 入 乳 熟 
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期 ,7 月 上 旬 开 始 收割 ,影像 NDVI 值 减 小 ,所 以 基于 
5 月 26 日 和 7 月 10 日 影像 的 NDVI 值 可 区 分 育 获 
和 小麦 。 

剩余 待 分 地 物 类 型 包括 玉米 、 天 然 草地 和 林地 
等 ,研究 区 土壤 盐 碱 化 程度 严重 ,天 然 草 地 以 盐 爪 爪 
等 耐 盐 植物 为 主 , 和 葡萄、 枸杞 等 植被 种 植 间距 较 大 ， 
其 生物 量 和 植被 覆盖 度 较 小 ,8 月 中 旬 玉 米 处 于 灌 
浆 期 ,NDVI 值 远 高 于 天 然 草 地 ,因此 可 基于 8 月 16 
日 影像 的 NDVI 值 ,区 分 玉米 和 天 然 草 地 以 及 葡萄 
和 枸杞 等 。 据 政府 相关 部 门 统计 报告 显示 ,在 10 月 
上 旬 研 究 区 玉米 已 基本 完成 收获 ,所 以 玉米 地 块 10 
月 12 日 遥感 影像 的 NDVI 远 低 于 林地 ,可 以 据 此 区 
分 玉米 和 林地 ,实现 玉米 种 植 范围 的 提取 。 

结合 作物 生长 物候 特征 与 作物 NDVI 和 NDWI 
时 序 特征 曲线 ,经 多 次 调试 构建 决策 树 分 类 规则 如 
图 4 所 示 : 

结合 作物 时 间 序 列 的 曲线 特征 和 耕作 信息 确 
立 分 类 时 间 ,选择 分 类 影像 和 特征 参数 构建 决策 
树 分 类 规则 ,操作 过 程 简单 ,分 类 规则 明了 ,无 需 
对 决策 树 分 类 规则 进行 剪 枝 , 但 其 分 类 时 相 的 选 
择 和 分 类 阔 值 的 判定 具有 较 大 的 主观 性 , 且 绘制 
作物 特征 曲线 时 采用 均值 , 易 导 致 部 分 特征 值 高 
于 或 低 于 平均 值 的 作物 被 漏 分 错 分 。CART 决策 
树 是 在 1984 年 由 Breiman 等 3 提出 的 一 种 二 又 
树 分 类 方法 , 它 能 够 充分 利用 空间 辅助 信息 ,对 影 
像 的 多 个 特征 变量 进行 选择 性 组 合 , 从 而 实现 对 
未 知 数据 的 分 类 、 预 测 和 挖掘 。 研 究 发 现 利用 二 
次 差分 法 计算 波峰 个 数 提取 首 禾 效果 较 好 ,已 有 
研究 的 阔 值 也 能 很 好 地 将 植被 与 非 植 被 区 分 开 
来 ,因此 只 对 水 稳 \、 玉 米 、 小 麦 和 其 他 植被 的 分 类 
规则 进行 重 构 。 通 过 目 视 解 译 联 合 Google Earth 
选取 包括 天 然 草地 .芦苇 和 葡萄 等 其 他 植被 的 样 
本 点 ,分 为 训练 样本 和 验证 样本 。 基 于 CART 决 
策 树 的 分 类 思想 ,将 训练 样本 作为 目标 变量 导入 
Clementine 软件 , 用 时 序 参数 作为 测试 变量 进行 学 


遥感 影像 时 间 序 列 数据 集 
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图 4 曲线 特征 决策 树 分 类 流程 及 参数 


Fig.4 Classification process and parameters of decision 


tree based on curve features 


习 并 剪 枝 ,建立 CART 决策 树 分 类 规则 ( 表 4) 。 
对 比 2 种 不 同方 法 建立 的 决策 树 ,除了 区 分 玉 
米 和 其 他 植被 时 , CART 决策 树 系 统 自动 选取 了 3 
月 26 日 的 NDVI 特征 ,而 基于 时 间 序 列 曲线 特征 建 
立 的 决策 树 选取 了 8 月 16 日 的 NDVI 特征, 其余 分 
类 时 相 选 择 均 相 同 ,对 应 的 分 类 阔 值 差 范围 仅 为 
0.056 ~0. 133 ,说 明 结 合 不 同 植被 的 生育 特点 ,依据 
不 同 作物 时 间 序 列 曲线 特征 选取 分 类 时 相 和 确定 分 
类 阔 值 是 可 行 的 。 绿 洲 作 为 干旱 区 具有 了 明显 小 气候 
效应 的 生态 景观 ,不 同 作 物 在 生育 期 内 物候 特征 差 
异 较 大 ,基于 遥感 提取 其 农作物 种 植 结 构 具 有 较 强 


表 4 CART 决策 树 分 类 规则 
Tab.4 Classification rules of CART decision tree 


条 件 结果 
NDVIw <0. 367 && NDVlogi >0. 534 Ih # 
NDWloss > -0.244 && NDWIio0 <0. 182 && NDVIomi >0. 367 KE 
NDVIw >0.367 && NDWIosz < -0.244 && NDVLon <0. 235 && NDWIowe <0. 177 玉米 


NDVI616 <0. 534 && NDW1 355 >0.177 && NDVI io >0.235 && NDWI, 9). >0. 182 其 他 
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干旱 区 研究 
的 典型 性 和 代表 性 。 麦 的 分 类 闭 值 设置 过 大 ,使 部 分 履 盖 度 较 高 的 天 然 
草地 被 误 分 为 小 麦 。 
4 结果 与 分 析 不 同 决策 树 方法 识别 的 不 同 作物 各 项 识别 精度 
均 较 高 。 首 蒂 的 提取 基于 时 序曲 线 波峰 个 数 的 计 
4.1 分 类 结果 


利用 ArcMap 对 提取 的 数据 进行 制图 ,得 到 研 
究 区 内 不 同 作物 的 空间 分 布 ( 图 5)。 
4.2 分 类 精度 检验 

基于 验证 样本 集 , 计 算 不 同 分 类 方法 的 混淆 和 矩 
EE ,对 研究 区 各 作物 分 类 结果 进行 验证 ( 表 5) 。 

MRS 中 可 以 看 出 ,2 种 决策 树 分 类 结果 都 较 
UE, CART 决策 树 总 体 分 类 精度 为 93. 06 ,Kappa 系 
数 为 0. 90 ,而 利用 曲线 特征 构建 的 决策 树 总 体 分 类 
精度 为 91.45 ,Kappa 系数 为 0.88。CART 决策 树 的 
总 体 分 类 精度 仅 比 基于 曲线 特征 构建 的 决策 树 高 
1.61 ,Kappa 系数 高 0.2。 对 比 不 同 分 类 方法 的 分 类 
结果 图 和 分 类 精度 可 发 现 , 基 于 时 间 序 列 曲 线 特征 
构建 的 决策 树 的 分 类 结果 中 少 部 分 玉米 被 误 分 为 水 
稳 ,小 麦 的 分 布 范围 和 面积 被 放大 ,主要 是 因为 在 依 
据 曲 线 特 征 确立 水 稻 分 类 规则 时 ,5 月 26 日 NDWI 
分 类 闭 值 偏 大 导致 部 分 玉米 被 错 分 为 水 稻 , 确立 小 
麦 的 分 类 规则 时 ,基于 曲线 特征 构建 的 决策 树 在 小 


算 , 其 制图 精度 和 用 户 精 度 均 达 到 95% 左右 ,说 明 
基于 二 次 差分 法 提取 曲线 极 值 进而 计算 波峰 个 数 识 
别 首 蒂 效 果 较 好 。 水 稻 在 不 同方 法 下 其 制图 精度 和 
用 户 精 度 均 在 91% 以 上 ,其 原因 主要 是 NDW 较 好 
地 区 分 出 了 水 稻 和 其 他 作物 , 另 一 方面 可 能 在 于 水 
稳 种 植 范围 连 片 集中 , 少 有 其 他 小 面积 作物 会 被 错 
分 到 水 稳 中 。 对 比 其 他 作物 ,玉米 的 分 类 精度 尤其 
是 用 户 精度 较 低 ,这 可 能 是 因为 玉米 的 播种 时 间 存 
在 差异 ,长 势 不 尽 相同 ,病人 值 范围 较 大 ,导致 与 玉米 
光谱 特征 相似 的 其 他 地 物 被 错 分 为 玉米 。 

此 外 ,影响 分 类 精度 的 原因 主要 体现 在 两 个 方 
面 ,首先 ,为 了 避免 平滑 处 理 影响 时 序曲 线 的 波峰 个 
数 提取 ,未 进行 时 间 序 列 数据 的 重 构 和 平滑 处 理 , 导 
致 未 能 剔除 遥感 数据 在 采集 过 程 中 受到 水 汽 、 气 溶 
胶 、 双 向 反射 及 数据 传输 等 多 重 因素 的 影响 而 出 现 
异常 植被 指数 ,在 一 定 程度 上 影响 了 提取 精度 。 其 
次 ,对 于 不 同 作物 交界 处 可 能 形成 的 混合 像 元 未 作 
进一步 的 分 类 ,也 会 使 分 类 精度 降低 。 
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图 $ 作物 种 植 结构 提取 效果 
Fig.5 Extracted results of crop planting structure 
表 5 不 同方 法 分 类 精度 比较 
Tab.5 Compared results of the classification accuracy of different methods 

Se 总 体 Kappa 制图 精度 用 户 精 度 
“精度 系数 水 稻 玉米 ME He 其 他 KE EK 小 麦 BRO 其 他 
线 特征 91.45 0.88 91.06 89.24 89.73 95.32 96.49 93.14 86.78 94.24 93.14 92.44 
CART 93.06 0.90 91.06 93.99 92.91 95.32 96. 49 96.26 87.19 94.93 95. 88 90.91 
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5 结论 


宁夏 河套 绿洲 地 势 较为 平坦 ,整个 生长 期 内 不 
同 作物 物候 特征 差异 较 大 , 且 在 提取 作物 类 型 之 前 
已 剔 除 水 体 和 其 他 建筑 物 .荒地 等 非 植被 的 影响 , 整 
体 分 类 结果 较为 理想 。 结 合 实地 调查 和 查阅 相关 统 
计 , 分 类 结果 与 实际 种 植 结构 吻合 。 

本 研究 利用 GF -1AWFV 数据 多 时 相 和 多 光谱 
EE ,构建 NDVI 和 NDWI 时 间 序 列 ,利用 不 同方 法 
构建 决策 树 ,实现 了 对 宁夏 河套 绿洲 作物 种 植 结构 
的 快速 .有效 .经济 的 遥感 提取 ,主要 结论 如 下 : 

(1) 不 同 作物 生育 期 不 同 ,时 间 序 列 曲线 形态 
也 不 同 ,基于 GF - 1/WEV 获得 的 归 一 化 植被 指数 
(NDVI) 时 间 序列 数据 。 能 够 描述 研究 区 农作物 不 
同时 期 的 生长 状态 ,表现 了 作物 在 季 相 节律 中 的 变 
化 过 程 ;基于 GF - 1/WEV 获得 的 归 一 化 水 体 指数 
NDWI 时 间 序 列 数据 ,能 够 描述 研究 区 农作物 不 同 
时 期 的 土壤 水 分 状态 ,表明 不 同 作物 灌溉 制度 的 差 
异 ;NDVI 时 间 序 列 及 NDWI 时 间 序 列 能 够 准确 表 
达 各 地 物 的 地 表 动 态 变化 信息 ,依据 作物 生长 周期 
中 的 关键 节点 选择 分 类 影像 及 特征 参数 ,结合 作物 
时 间 序 列 曲线 特征 ,可 以 有 效 地 识别 不 同 作 物 。 

(2) 2 种 决策 树 分 类 方法 分 类 精度 均 较 高 ,但 
基于 曲线 特征 构建 决策 树 操 作 过 程 更 简单 ,规则 更 
明了 , 且 无 需 对 决策 树 进 行 剪 枝 , 而 CART 决策 树 能 
人 够 更 精确 地 确立 各 作物 分 类 阔 值 ,有 效 提 高 分 类 
精度 。 

(3) 对 比 CART 决策 树 ,结合 不 同 植被 生育 特 
点 ,基于 主要 作物 曲线 特征 构建 的 决策 树 在 分 类 时 
相 、 分 类 阔 值 和 分 类 精度 都 相差 不 大 ,说明 研究 区 内 
不 同 作物 在 生育 期 内 物候 特征 差异 较 大 ,其 对 应 的 
NDVI 和 NDWI 时 间 序 列 曲线 特征 差异 明显 ,对 于 
具有 明显 小 气候 效应 的 干旱 区 绿洲 ,利用 时 间 分 辩 
率 和 空间 分 辨 率 都 较 优 的 GF -1 - WEV 时 间 序 列 
数据 ,对 其 农作物 种 植 结构 进行 愧 感 提取 具有 较 强 
的 实用 性 和 代表 性 。 
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Extraction of Crops in Oasis Based on GF -1/WFV Time Series 


LIU Ya-qing'*, WANG Lei'?”, ZHAO Xi-ni'?, QU Xiang-ning'*, XU Xing, WANG Rui!” 
(1. Breeding Base for State Key Laboratory of Land Degradation and Ecosystem Restoration in Northwest China, Ningxia University , 
Yinchuan 750021 , Ningxia , China; 

2. Key Laboratory for Restoration and Reconstruction of Degenerated Ecosystem in Northwest China under Ministry of Education , 
Ningxia University , Yinchuan 750021 , Ningxia , China ; 

3. International Institute for Earth System Science , Nanjing University , Nanjing 210093 , Jiangsu , China) 


Abstract: The rapid extraction of regional crops is of great significance for agricultural production management, 
planting structure adjustment and optimization, and food security. The use of time series data of remote sensing im- 
age for extracting crops is an important means, and the time resolution and spatial resolution are constraints. The 
time resolution of medium spatial resolution of remote sensing data is significantly improved with the successful 
launch of the first satellite GF — 1 of China High-resolution Earth Observation System. However, the seasonal rhythm 
features of crops differ greatly under different climate backgrounds. The Hetao Oasis is irrigated by the Yellow Riv- 
er. It has an ecological landscape with obvious microclimate effect. The extraction of remote sensing data of crops is 
typical and representative. The Hetao Oasis in Ningxia is selected as the study area to analyze and evaluate the ap- 
plicability of GF -1 satellite data in extracting crops. According to the characteristics of crop planting in the study 
area, the four main crops including rice ,corn, wheat and alfalfa were extracted. Firstly ,the GF -1/WFV image was 
calculated by band,and the NDVI and NDWI time series of remote sensing data sets were obtained. The values of 
NDVI and NDWI of the main crops were extracted, and the time series curves of NDVI and NDWI were construc- 
ted. The time series data of remote sensing indexes of the main crops were analyze ,the farming information was de- 
termined to extract the phases and thresholds of the crops, construct the decision tree and extract the crops. The time 
series data of natural grasslands,reed and grape and of other vegetation types were combined to establish the CART 
decision tree so as to verify the classification effect of the correlation time series and the decision tree constructed by 
the farming information. The results are as follows : Q) The morphological characteristics of time series curves were 
different from different farming ways,and the NDVI time series data could be used to describe the growth status of 
crops in different periods; the NDWI time series data obtained from GF -1/WFYV could be used to describe the soil 
moisture status of different crops and reflect the different irrigation systems of different crops in the study area. The 
NDVI time series and the NDWI time series could be used to accurately express the surface dynamic change infor- 
mation of the surrounding areas of the study area, correlate the crop time series with farming information, and effec- 
tively recognize different crops;(2) Compared with the CART decision tree , the decision tree of classification phase , 
classification threshold and classification accuracy constructed by combining different crop farming information and 
main crop curve features were similar. These indicated that the farming ways were different from different crops in 
the study area, and the time series of corresponding NDVI and NDWI values were significantly different. 

Key words: crops; cropping pattern; GF -1/WFYV time series data; decision tree; remote sensing mapping; 


oasis; Ningxia 
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